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因 累 积分 二 次 项 与 母亲 消极 教养 的 交互 
影响 两 个 时 间 点 的 青少年 攻击 行为 : 在 较 高 母亲 消极 教养 条 件 下 ,携带 较 多 或 较 少 低 多 巴 胺 
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条 件 下 ,多 基因 累积 分 二 次 项 与 青少年 的 攻击 行为 关系 不 显著 。 本 研究 为 多 巴 胺 系统 基因 的 
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1 引言 

攻击 是 个 体 有 意 对 他 人 身体 或 心理 
24.3% 的 中 国 青少年 具有 较为 严重 的 攻击 或 暴力 和 
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的 普 ; 


机 制 ， 淮 试 提出 个 体 化 的 攻击 行为 家 庭 预 防 方案 。 


年 攻击 行为 不 仅 导致 
大 的 社会 经 济 和 精神 
性 与 危害 性 , 研究 者 从 遗传 因素 和 家 庭 教 养 广 


青 少 


进行 伤害 的 行为 (Crick & Grotpeter, 1995) , 7.48% ~ 
7 为 (Wang et al., 2012; Yu et al., 2021) 。 

串 精 神 病 风险 和 暴力 犯罪 行为 , 还 会 造成 巨 
卫生 医疗 负担 (Ansary et al., 2017; Lopez et al., 2018) 。 由 于 攻击 行为 


1.1 多 巴 胺 系统 基因 与 攻击 
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用 探索 攻击 行为 发 生 的 基因 一 环境 作用 


因 参 与 情绪 调控 、 执行 控 和 


出 以 及 动机 过 程 (Gatzke-Kopp, 2011; Nikolova et 


al., 2011; Tielbeek et al., 2018; Vrantsidis et al., 2021) ， 与 个 体 的 攻击 行为 密切 相关 (Albaugh 


et al., 2010; van Goozen et al., 2016; Zhang et al., 2016) 。 然 而 ， 现 有 关于 多 巴 胺 系统 基因 与 
为 间 关 联 的 分 子 遗 传 研究 结果 存在 不 一 致 。 部 分 研究 显示 , 与 低 多 巴 胺 活性 相关 基因 
加 相关 联 (Hygen et al., 2015; Jones et al., 2001; van Goozen et al., 2016) 。 
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严重 4 


KAS Bad 
FE 活 事件 时 表现 更 多 攻击 行为 。 然 而 ， 也 有 研究 显示 Met 等 位 基因 与 高 攻击 行为 密切 


HA 


J 为 有 关 ，Hygen ÆA (2015) 的 研究 发 现 该 位 点 Val 纯 合 子 携带 者 经 历 较 多 
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击 行为 存在 关联 ; Zhang 等 人 (2016) 发 现 Met 等 位 基因 携带 者 在 经 历 低 水 平 积极 教养 时 表 
现 出 更 多 反应 性 攻击 。 尽管 研 究 者 试图 从 样本 特征 、 环境 因素 以 及 攻击 测评 方式 角度 对 研究 
结果 的 分 歧 进行 解释 ， 但 是 既 有 研究 可 能 忽视 了 两 个 重要 事实 : (1) 多 巴 胺 系统 单个 基因 
作用 微 效 , 多 个 基因 存在 联合 与 交互 作用 , 多巴胺 系统 遗传 效应 积分 与 攻击 行为 并 非 线性 关 
A: (2) 多 巴 胺 系统 基因 的 作用 受 环 境 因 素 的 调节 。 

神经 化 学 和 心理 遗传 研究 表明 多 巴 胺 系统 基因 与 攻击 行为 间 可 能 存在 非 线性 关系 ,神经 
生化 研究 证 据 表 明 ， 攻 击 行为 相关 脑 功能 与 多 巴 胺 活性 间 呈 倒 U 型 关系 ， 中 等 多 巴 胺 活性 
水 平 下 PFC 及 相关 脑 区 的 功能 最 优 ， 多 巴 胺 活性 过 高 或 过 低 时 PFC 及 相关 脑 区 功能 受 损 
(Arnsten, 2009; Bertolino et al., 2009; Robbins & Arnsten, 2009; Seamans & Yang, 2004; Tian 
et al., 2013; Tunbridge et al., 2006) 。 心 理 遗 传 研究 发 现 多 巴 胺 系统 基因 对 攻击 行为 相关 脑 区 
结构 和 功能 的 影响 存在 交互 作用 或 非 线性 关系 。 例 如 ，Bertolino A (2009) 探究 了 DAT, 
DRD2 基因 与 特定 脑 区 活性 的 关联 ， 发 现 携带 与 高 多 巴 胺 活性 相关 的 基因 型 组 合 (DAT 
10R/10R 和 DRD2 GG 基因 型 ) 或 携带 与 低 多 巴 胺 活性 相关 的 基因 型 组 合 (D4T9R 和 DRD2 
GT 基因 型 ) 个 体 的 尾 状 核 和 额 中 回 活 性 较 低 、 灰 质 体积 更 大 ， 而 携带 中 等 多 巴 胺 活性 相关 
基因 型 组 合 个 体 的 尾 状 核 和 额 中 回 活 性 、 灰 质 体积 处 于 正常 水 平 。 因 此 ， 我 们 预期 ， 多 巴 胺 
系统 基因 与 攻击 行为 间 可 能 存在 非 线 性 关系 。 
1.2 多 巴 胺 系统 多 基因 的 联合 效应 

近年 来 ， 基 于 多 基因 假说 的 多 基因 累积 分 (multilocus genetic profile score, MGPS 
受到 重视 。 该 方法 通过 累计 多 个 基因 位 点 获得 多 基因 累积 分 , 揭示 多 个 基因 对 某 种 心理 或 行 
为 个 体 差 异 的 贡献 (Cao & Rijlaarsdam, 2022; Davies et al., 2019; Nikolova et al., 2011; 
Thibodeau et al., 2015) 。 

神经 系统 内 多 巴 胺 活性 受到 多 巴 胺 能 水 平 、 受 体 效 能 和 利用 效率 的 影响 。 因 此 ， 参 与 多 
巴 胺 降解 、 传 导 与 转运 的 基因 间 的 联合 效应 共同 调节 突 触 间隙 内 的 多 巴 胺 活性 , 影响 儿童 青 
少年 心理 病理 问题 (Davies et al., 2019; Nikolova et al., 2011; Thibodeau et al., 2015) 。 儿 茶 酚 
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胺 氧 位 甲 基 转 移 酶 〈catechol-O-methyltransferase, COMT) ~ D2 型 多 巴 胺 受 体 (D2 receptor, 
DRD2) 和 多 巴 胺 转运 体 (dopamine transporter, DAT1) 基因 分 别 影响 多 巴 胺 的 降解 、 传 导 
与 转运 (Chester & DeWall, 2016; Nikolova et al., 2011) ， 共 同调 节 多 巴 胺 活性 。COMT 基因 
编码 儿 茶 酚 胺 氧 位 甲 基 转 移 酶 , 负责 降解 多 巴 胺 等 单 胺 类 神经 递 质 。 COMT 184680 多 态 位 点 
翻译 产物 中 第 158 (AER REST AIR (valine, Val) 到 和 蛋氨酸 (methionine, Met) 的 替换 使 
COMT 酶 活性 降低 40% (Chen et al., 2004; Lotta et al., 1995; Meyer et al., 2016) . DRD2 基因 
负责 编码 D2 型 多 巴 胺 受 体 。 该 基因 启动 子 区 rs1799978 (A-241G) 多 态 性 在 -241 位 置 处 发 
ERRI CA) ASEE (G) 的 基因 转换 ， 与 更 多 的 多 巴 胺 D2 受 体 密度 和 转录 水 平 有 关 ， 
进而 影响 多 巴 胺 活性 (Genis-Mendoza et al., 2018; Xing et al., 2007; Zhang & Malhotra, 2011). 
DATI 基因 负责 编码 多 巴 胺 转运 体 。 该 基因 3? 端 非 翻 译 区 rs27072 位 点 工 等 位 基因 比 C 等 位 
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基因 型 具有 更 高 转 率 水 平 ， 从 而 在 突 触 中 更 有 效 地 重新 摄取 多 巴 胺 , 减少 突 触 间隙 的 多 巴 胺 
JKF (Pinsonneault et al., 2011) 。 这 三 种 基因 在 前 额 叶 、 纹 状 体 和 中 脑 等 攻击 相关 脑 区 中 表 
达 水 平 更 高 (Yamaguchi & Lin, 2018) 。 已 有 研究 发 现 ， 这 三 个 基因 多 态 性 位 点 与 儿童 青 少 
年 攻击 行为 密切 相关 (Chester & DeWall, 2016; Davies et al., 2015; van Goozen et al., 2016; Zai 


et al., 2012; Zhang et al., 2016), 协同 影响 攻击 行为 (Davies et al., 2019; Vrantsidis et al., 2021). 
本 研究 拟 选 取 COMT, DRD2 和 DATI 的 上 述 三 种 高 功能 等 位 基因 创建 多 基因 累积 分 ， 考 察 
多 基因 累积 分 与 攻击 行为 的 关联 。 
1.3 母亲 消极 教养 的 调节 作用 

多 巴 胺 系统 基因 对 攻击 行为 的 影响 作用 受 环境 因素 调节 ， 即 存在 G x E 交互 作用 
(Beauchaine et al., 2007; Davies et al., 2019; Shanahan & Hofer, 2005) 。 生 物 发 展 模型 
(biosocial developmental model, BDM) 指出 ， 重 复 暴 露 于 不 利 教养 环境 会 增加 携带 多 巴 胺 
功能 缺陷 基因 个 体 发 展 出 问题 行为 的 风险 (Beauchaine et al., 2007) 。 近 期 ，Tomlinson 等 
(2022) 的 双生 子 研究 为 此 提供 了 证 据 支持 ， 他 们 发 现 ， 温情 (积极 教养 ) 减弱 遗传 因素 对 
冷酷 无 情 特质 的 影响 ， 但 是 遗传 素质 对 冷酷 无 情 的 影响 在 严厉 管教 〈 消 极 教养 ) 中 更 强 。 

家 庭 教养 对 儿童 青少年 攻击 行为 具有 重要 影响 (Davies et al., 2019; Gershoff, 2002; 
Stocker et al., 2017; Thibodeau et al., 2015) 。 其 
教养 行为 会 使 得 儿童 青少年 出 现 对 不 良 行为 的 模仿 、 增 强 敌 意 情绪 情感 及 认 知 归 因 等 
(Gershoff, 2002; Goulter et al., 2020) ， 限 制 、 阻 碍 情绪 及 行为 调节 功能 发 展 等 〈Ganiban et 
al., 2021; Thibodeau et al., 2015) ， 进 而 促进 攻击 行为 发 生 (Zhang et al., 2016) 。 既 有 研究 
揭示 ， 消 极 教养 与 多 巴 胺 系统 多 基因 交互 影响 儿童 青少年 攻击 相关 行为 ， 如 反 社 会 行为 、 外 
化 问题 等 (Davies et al., 2015, 2019; Thibodeau et al., 2015) 。 璧 如 ，Davies 等 (2015, 2019) 
的 研究 发 现 ， 多 巴 胺 系统 多 基因 累积 分 与 母亲 非 敏感 性 等 消极 教养 交互 预测 个 体 的 外 化 问 
题 ,携带 与 较 低 多 巴 胺 活性 相关 等 位 基因 的 个 体 在 经 历 高 水 平 的 消极 教养 时 表现 出 相对 较 高 
水 平 外 化 问题 , 经 历 低 水 平 消极 教养 时 外 化 问题 水 平 相对 较 低 。 因 此 ， 本 研究 中 ， 为 揭示 多 
巴 胺 系统 基因 与 攻击 行为 的 关联 , 拟 选择 母亲 消极 教养 作为 环境 因素 , 考察 其 与 多 巴 胺 系统 
多 基因 对 攻击 行为 的 交互 作用 , 通过 揭示 高 消极 教养 条 件 下 多 巴 胺 系统 多 基因 与 攻击 行为 间 
的 关联 阐明 多 巴 胺 系统 基因 对 攻击 行为 的 影响 。 
1.4 本 研究 问题 与 假设 

青少年 期 (adolescence) 是 个 体 发 展 中 情绪 情感 、 认 知 以 及 神经 生物 功能 发 生 重要 变化 
的 时 期 。 青 少年 期 PFC 内 基线 多 巴 胺 含量 增加 (Naneix et al., 2012) ， 脑 内 多 巴 胺 受 体 表达 
(Naneix et al., 2012) 以 及 多 巴 胺 转运 体 密度 (Tielbeek et al., 2018) 逐渐 增加 并 达到 峰值 ， 
因而 ， 青 少年 期 多 巴 胺 能 敏感 性 增强 ， 影 响 多 巴 胺 调节 通路 功能 。 青 少年 期 多 巴 胺 通路 调节 
功能 的 改变 使 PFC、 纹 状 体 的 结构 与 功能 ， 以 及 PFC 与 杏仁 核 间 功 能 连接 发 生 显著 变化 
(Romeo, 2017; Toth, 2022; Tottenham & Galvan, 2016) ， 个 体 表现 出 对 奖惩、 威胁 性 刺激 等 


敏感 性 与 反应 性 增强 (Tielbeek et al., 2018) 。 由 于 青少年 期 多 巴 胺 能 敏感 性 增强 ， 个 体 所 
携带 的 与 高 或 低 多 巴 胺 活性 相关 等 位 基因 的 作用 可 能 会 被 进一步 增强 , 使 多 巴 胺 系统 基因 和 与 
行为 间 关联 表现 出 不 同 于 其 他 发 展 阶段 的 模式 。 

既 有 关于 基因 及 其 与 环境 交互 对 攻击 行为 影响 的 研究 主要 聚焦 于 幼儿 或 童年 中 晚期 
(Davies et al., 2015; Thibodeau et al., 2015; Vrantsidis et al., 2021) ， 少 有 研究 探讨 青少年 期 
(Davies et al., 2019) 。 本 研究 将 基于 多 基因 研究 视角 ， 以 COMT. DATI 和 DRD2 三 个 基 

因 多 态 性 累积 分 为 遗传 指标 , 母亲 消极 教养 为 环境 指标 , 青少年 期 2 次 测量 的 攻击 作为 攻 
行为 指标 ， 探 究 多 巴 胺 系统 多 基因 与 青少年 期 攻击 行为 的 关系 以 及 母亲 消极 教养 的 调节 作 
Jo 基于 已 有 研究 , 我 们 拟 探索 多 基因 累积 分 与 青少年 攻击 行为 间 的 非 线 性 关系 以 及 母亲 教 
养 在 其 中 的 调节 作用 ,我们 预期 携带 与 高 或 低 多 巴 胺 活性 相关 等 位 基因 的 青少年 均 可 能 表现 
出 高 水 平 攻 击 行为 , 母亲 消极 教养 调节 多 基因 累积 分 与 青少年 攻击 行为 间 的 关系 , 母亲 消极 
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= 教养 水 平 较 高 时 多 基因 累积 分 与 青少年 攻击 行为 间 的 非 线性 关系 更 突出。 
© 2 研究 方法 
一 2.1 被 二 


本 研究 被 试 来 自 一 项 大 型 追踪 项 目 ， 是 利用 该 追踪 项 目 被 试 开展 的 遗传 研究 。 本 研究 
中 ， 初 一 学 年 末 《T1) 测量 母亲 消极 教养 行为 ， 初 二 学 年 末 〈T2) 收集 基因 数据 ，T1、T2 
两 次 收集 青少年 攻击 行为 数据 。T1、T2 为 期 1 年 的 时 间 间 隔 遵循 纵向 研究 设计 常规 模式 ， 
该 时 间 间 隔 能 反映 出 有 意义 的 发 展 变化 或 差异 性 , 不 会 因 间隔 过 短 使 重复 测评 无 意义 或 因 间 
隔 过 长 而 漏 评 有 意义 的 发 展 变化 、 差 异性 或 预测 效应 〈Collins, 2006) 。 此 外 , 在 同一 年 〈 学 
年 ) 类 似 时 间 重 复 测 量 ， 可 避免 季节 、 开 学 时 长 、 学 年 中 不 同 周 期 性 事件 等 因素 对 结果 变量 
的 影响 (Dormann & Griffin, 2015) 。 

受 研 究 经 费 所 限 ， 该 追踪 项 目 仅 对 部 分 被 试 进行 唾液 采集 和 DNA 分 型 。 具 有 遗传 数据 
(N = 1078) 与 没有 遗传 数据 CN = 927) 的 被 试 在 主要 人 口 学 信息 上 无 显著 差异 〈 年 龄 : 
人 2003)1 =—0.74, p = 0.46; SES?: (1867) = —0.95, p = 0.34; 性 别 : 2 =3.39, df= 1,p=0.07)， 但 
是 具有 遗传 数据 被 试 的 攻击 行为 与 消极 教养 较 低 〈T1 攻击 : (975) = 7.45, p < 0.001; T2 攻击 : 
(1042) = 6.24, p < 0.001; 消极 教养 : 1880) = 2.07, p= 0.04) 。 上 有 具有 遗传 数据 被 试 含 汉族 CN 
= 1044, 96.8%) 和 其 他 民族 CN = 34, 3.2%) 。 昌 然 汉 族 与 其 他 民族 被 试 在 基本 研究 变量 
性别、 基因 型 分 布 、 年 龄 、SES、T1 攻击 行为 、T2 攻击 行为 、 母 杀 消 极 教养 ) 无 显著 差 
异 Gs < 4.45, ps > 0.49; ltls < 1.58, ps > 0.12) ， 但 考虑 到 基因 功能 和 连锁 不 平衡 〈linkage 
diequilibrium, LD) 对 适应 结果 的 影响 存在 种 族 差 异 〈Avinun et al., 2020) ， 为 避免 由 于 民族 


ant 


pey 


1 独立 样本 1 检验 中 ,年龄 、SES、 攻击 行为 及 消极 教养 等 方差 不 齐 性 ,为 校正 统计 偏差 , 采用 校正 后 的 自由 度 计算 t 值 (Tabachnick 
& Fidell, 2019) 。 此 处 报告 度 是 校正 后 的 自由 度 。 

2 SES 指标 由 父母 职业 、 文 化 程度 及 家 庭 收入 合成 。 参 照 既 有 研究 (Chen etal., 2019) ， 父 母 职 业 根据 职业 的 专业 技术 性 程度 归 
为 三 类 : 农民 或 下 岗 失 业 人 员 、 蓝 领 、 专 业 或 半 专 业 性 人 员 ， 分 别 赋值 1 ~ 3。 父 母 受 教育 水 平 包括 初中 ( 含 未 毕业 ) 及 以 下 、 高 
中 或 中 专 ( 含 未 毕业 ) 、 大 学 〈 含 专科 ) 及 以 上 三 类 ， 分 别 赋值 1 ~3。 家 庭 月 收入 为 3000 元 以 下 、3000 元 到 6000 元 、6000 元 以 
上 ， 分 别 赋值 1 ~3。 参 照 已 有 研究 ( 范 兴 华 等 , 2012) ， 采 用 因子 分 析 法 合成 SES。SES 得 分 越 高 ， 代 表 家 庭 社会 经 济 地 位 越 高 。 


差异 造成 偏差 并 达到 最 大 统计 检验 力 (Cardon & Palmer, 2003) ， 参 照 已 有 研究 (如 , Avinun 
et al., 2020; Cao & Rijlaarsdam, 2022; Stocker et al., 2017) ， 本 研究 仅 对 汉族 被 试 进行 分 析 。 

最 终 样本 1044 人 ，T1 平均 年 龄 为 13.32 + 0.49 岁 ， 男 生 520 A (49.8%) ， 女 生 524 
A (50.2%) 。 母 亲 和 父 杀 受 教育 水 平 在 初中 及 以 下 者 分 别 占 11.6% 和 8.3%， 高 中 或 中 专 者 
分 别 占 26.5% 和 20.6%， 大 学 及 以 上 者 分 别 占 61.9% 和 71.1%。 家 庭 月 收入 低 于 3000 元 的 占 
23.1%，3000 ~ 6000 元 的 占 48.0%， 高 于 6000 元 的 占 28.9%。 
2.2 研究 程序 与 研究 工具 

研究 通过 本 单位 伦理 委员 会 审核 批准 。 数 据 收集 前 ， 将 问卷 施 测 、 唾 液 采集 流程 等 告知 
青少年 所 在 学 校 、 监 护 人 及 青少年 本 人 ， 取 得 知情 同意 。 

青少年 攻击 行为 数据 采用 问卷 法 ， 以 班级 为 单位 集体 施 测 ,结束 后 当场 收回 。 青 少年 中 
液 样本 以 班级 为 单位 组 织 采 集 。 唾 液 样 本 采集 严格 按照 规范 进行 , 采样 时 被 试 无 发 烧 ， 采 和 集 
前 30 分 钟 无 进食 、 吸 烟 、 饮 酒 、 饮 水 或 嚼 口香糖 ， 采 集 完 成 后 检查 以 确保 样本 质量 ， 置 于 
-20 °C 低温 保存 ， 一 周 内 冷 链 运输 至 生物 公司 提取 分 型 。 攻 击 行为 问卷 施 测 、 唾 液 样 本 采集 
时 , 每 班 由 2 名 经 过 严格 培训 的 发 展 心理 学 教师 和 研究 生 任 主 试 。 母亲 消极 教养 数据 由 青 少 
年 将 装 有 母亲 教养 问卷 的 信封 带 回 家 由 母亲 填写 ， 第 二 天 带 回 交 班 主任 老师 统一 收回 。 
2.2.1 多 巴 胺 系统 基因 提取 、 分 型 与 计 分 

青少年 唾液 样本 每 人 不 少 于 2 ml。 利 用 Sequenom (San Diego, CA, USA) 世 片 基质 加 
助 激光 解吸 /电离 飞行 时 间 (MALDI-TOF) 质谱 平台 ， 按 照 标准 化 技术 规范 、 流 程 对 唾液 样 
本 进行 DNA 提取 与 分 型 (Cao & Rijlaarsdam, 2022; Zhang et al., 2016) 。 
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COMT Val158Met(y? = 2.06, p = 0.15; Val/Val 56.8%, Val/Met 36.0%, Met/Met 7.2%) ,DRD2 
1s1799978 (y? = 0.004, p = 0.95; AA 66.5%, AG 30.1%, GG 3.4% ) 和 DATI rs27072 Q= 1.63, p= 
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= 0.20; TT 5.9%, CT 33.9%, TT 60.2%) 三 个 基因 多 态 性 的 基因 型 分 布 均 符合 Hardy-Weinberg 平 


衡 。 多 基因 研究 基于 “常见 疾病 一 常见 基因 (common disease-common variant) ”假设 ， 要 
求 次 要 等 位 基因 频率 (minor allele frequency, MAF) 大 于 5% (Manolio et al., 2009) . Att 
究 中 ， 三 个 基因 位 点 的 次 要 等 位 基因 频率 均 符合 要 求 《COMT: 25.2%, DRD2: 18.5%, DATI: 
22.9%) 。 三 个 多 态 性 的 基因 型 分 布 均 不 存在 显著 性 别 差 异 GPs < 1.24, df=2,ps>0.54) 。 
为 确定 COMT、DRD2 和 DAT] 三 个 基因 是 否 适 合 线性 累加 ， 进 行 系列 统计 检验 。 首 先 ， 由 
于 基因 编码 方式 会 影响 多 基因 累积 分 结果 ， 为 检验 线性 编码 方式 〈 即 0，1，2) 是 否 适用 于 形 
成 多 基因 累积 分 ， 参 照 已 有 研究 (Stocker et al, 2017; HHT, 张 文 新 , 2019) ， 进 行 线性 基因 
模型 检验 。 结 果 发 现 ， 对 T1、T2 攻击 ， 虽 然 分 解 模型 估计 的 参数 最 多 ， 但 其 解释 率 并 未 显著 高 
于 线性 基因 模型 〈 即 限定 同一 基因 的 显 性 和 隐 性 编码 估计 系数 相等 ; AR? = 0.004, 0.007, ps > 
0.28) ， 因 而 三 个 基因 位 点 的 各 基因 型 符合 线性 编码 方式 。 其 次 ， 为 了 排除 单个 基因 的 主导 效 
悄 ， 进 行 等 基因 效应 模型 检验 。 结 果 显 示 ， 等 基因 效应 模型 与 线性 基因 效应 模型 间 没 有 显著 
差异 〈AR2 = 0.002、0.007,ps > 0.11) ， 表 明 三 个 基因 对 攻击 行为 的 预测 强度 无 差异 ， 不 存在 


单个 基因 预测 力 过 强 主导 多 基因 累积 分 的 现象 。 因 此, 根据 个 体 携带 与 低 多 巴 胺 活性 相关 的 
等 位 基因 (COMT: Val/Val = 0, Val/Met = 1, Met/Met = 2; DAT: AA = 0, AG = 1, GG = 2; DRD2: 
TT = 0, CT= 1, CC =2) 的 数量 进行 线性 编码 ， 通 过 求 和 获得 基因 位 点 累积 分 。 多 巴 胺 系统 
多 基因 累积 分 数 分 布 为 0 (N= 24, 2.3%) 、1 (N= 164, 15.7%) 、2 (N= 378, 36.2%) 、3 
(N=334, 32.0%) . 4 (N=120, 11.5%) 和 5 (N=24,2.3%) 。 
2.2.2 母亲 消极 教养 行为 

采用 中 文 版 儿童 养育 实践 问卷 (Child-Rearing Practices Report, CRPR; Chen et al., 2000), 
由 母 杀 报告 母亲 消极 教养 行为 。 该 问卷 在 中 国 儿 童 青少年 中 广泛 使 用 , 具有 良好 的 信 效 度 ( 如 ， 
Zhang et al., 2017) 。 在 本 研究 中 ， 惩 罚 维度 的 1 个 条 目 “ 我 告诉 我 的 孩子 ， 如 果 他 〈 她 ) 
做 了 坏事 ， 就 要 受到 惩罚 ”与 消极 教养 总 分 相关 低 ， 参 照 已 有 研究 (Zhang et al, 2017) 将 
其 删除 。 本 研究 采用 拒绝 (4 条目， 如“ 如果 我 的 孩子 不 来 烦 我 ， 我 就 不 会 理 皮 她 ”) 、 惩 
Tl (6 条 目 ， 如 “我 认为 体罚 是 管教 孩子 的 最 好 的 方式 ”) 维度 共 10 条 目 均 分 作为 母亲 消 
极 教养 行为 指标 。 问 卷 采 用 5 点 计 分 〈 从 “完全 不 符合 ”= 0 到 “完全 符合 ”= 4) ， 得 分 越 
高 表明 母亲 消极 教养 行为 越 多 。 验 证 性 因素 分 析 显 示 ， 单 因子 模型 拟 合 良 好 GO? = 152.73, df 
= 35, RMSEA = 0.058, CFI = 0.91, TLI = 0.89)。 本 研究 中 母亲 消极 教养 行为 测量 的 Cronbach’s 
a 系数 为 0.72。 
2.2.3 攻击 行为 

采用 同伴 提名 法 测量 青少年 攻击 行为 (Crick & Grotpeter, 1995) 。 该 量 表 在 青少年 攻击 

究 中 广泛 应 用 ， 能 够 提供 攻击 行为 有 效 可 靠 的 数据 (Mehari etal., 2019) ， 具 有 和 良好 信 效 
度 。 问 卷 包含 身体 攻击 〈 如 “经 常 打 人 的 同学 ”) 和 关系 攻击 〈 如 “ 谁 经 常 在 背后 说 别 的 同 
学 的 坏话 ”) 方面 的 4 个 条 目 。 被 试 按照 每 个 题目 描述 的 行为 ， 从 班级 中 选 出 自己 认为 最 符 
合 描述 的 三 名 同学 。 合计 每 个 被 试 得 到 的 提名 频次 并 进行 班 内 标准 化 , 作为 被 试 攻击 行为 得 
分 。 两 个 时 间 点 攻击 行为 测量 的 Cronbach’s a 系数 分 别 为 0.88 和 0.87. 
2.3 数据 处 理 与 分 析 

首先 采用 相关 分 析 考 察 研 究 变量 间 的 关系 ， 其 次 采用 SPSS 26.0 进行 分 层 线性 回归 分 析 
考察 多 基因 累积 分 与 母亲 消极 教养 对 青少年 攻击 行为 的 影响 。 分 层 回 归 分 析 模 型 第 一 层 纳入 
性 别 作为 控制 变量 ， 第 二 层 纳入 多 基因 累积 分 、 母 亲 消 极 教养 为 预测 变量 ， 第 三 层 为 多 基因 
累积 分 与 母亲 消极 教养 的 交互 项 ， 第 四 层 、 第 五 层 分 别 纳入 多 基因 累积 分 二 次 项 、 其 与 母亲 
消极 教养 的 交互 项 。 使 用 包含 二 次 项 (本 研究 中 ,为 多 基因 累积 分 二 次 项 ) 的 回归 分 析 能 够 
揭示 常见 非 线性 关系 ， 含 U 型 、 倒 U 型 ， 以 及 预测 变量 对 结果 变量 的 预测 随 预 测 变 量 增 大 
而 增强 或 减弱 的 其 他 非 线性 关系 模式 ， 二 次 项 显著 且 为 正则 为 U 型 关系 ， 二 次 项 显著 且 为 
负 则 为 倒 U 型 关系 (Cohen et al., 2003; Haans et al., 2016) 。 若 任 一 交互 项 显著 ， 则 进行 简 
单 斜率 检验 (Cohen et al., 2003) 。 简 单 斜 率 检 验 的 做 法 是 ， 按 照 +1 SD 将 母亲 消极 教养 区 
分 为 高 组 和 低 组 , 分 析 高 、 低 母亲 消极 教养 条 件 下 多 基因 累积 分 、 二 次 项 与 攻击 行为 的 关联 ， 
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同时 根据 既 有 研究 建议 (Kochanska et al., 2011; Roisman et al., 2012; Zhang et al., 2015) 按照 
+2 SD 将 母亲 消极 教养 区 分 为 高 组 和 低 组 , 以 探索 和 阐明 表现 于 +1 SD 之 外 的 显著 交互 作用 。 
为 避免 多 重 共 线 性 ， 多 基因 累积 分 和 母亲 消极 教养 均 进行 标准 化 处 理 。 

为 验证 结果 的 可 靠 性 进行 系列 敏感 性 分 析 : (1) 对 任 一 单 基因 重复 上 述 回 归 模 型 检验 ， 
比较 单 基 因 与 多 基因 累积 分 模型 最 著 性 的 差异 (Cao & Rijlaarsdam, 2022); (2) 使 用 MATLAB 
R2013b 的 regress 和 crossvalind 函数 ， 进 行 k fold 交叉 验证 (Marcot & Hanea, 2021; sKIC# 
等 , 2021) ; G) 将 样本 按照 50% 随 机 分 为 两 个 子 样本 进行 上 述 回归 分 析 ， 以 进行 内 部 验 
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(Cao & Rijlaarsdam, 2022; Zhang et al., 2016) 。 


3.1 描述 统计 结果 
变量 描述 统计 与 相关 系数 见 表 1。 多 基因 累积 分 与 母亲 消极 教养 相关 不 显著 ， 排 除了 
基因 一 环境 相关 。 多 基因 累积 分 与 两 个 时 间 点 的 攻击 相关 不 显著 , 母亲 消极 教养 与 两 时 间 点 
的 攻击 行为 均 显 著 正 相关 。 独 立 样本 上 检验 表明 ， 两 个 时 间 点 的 攻击 行为 、 母 亲 消 极 教养 行 
为 性 别 差 异 显著 (T1 攻击 3: (782) = -7.13, T2 攻击 : (772) = -6.68， 消 极 教养 : {(958) = -3.65， 
ps<0.001) ， 男 生 的 攻击 行为 和 母亲 消极 教养 行为 显著 高 于 女生 。 

表 1 描述 统计 与 相关 系数 


ZE MM D ”性别 * 多 基因 累积 分 消极 教养 ”TI 攻击 T2 攻 击 ， 
性 别 一 = 1 
多 基因 累积 分 一 一 -0.02 1 
消极 教养 0.99 0.48 0.11 0.03 1 
Tl 攻击 —0.12 0.68 0.22*** -0.02 0.12*** 1 
T2 攻击 -0.09 0.71 0.20*** -0.01 0.12*** 0.76 1 


TE: MRETI, SDRE. EMA, O- 女生 ，1 BE. p< 0l, p<005 "p <00, p<0001, Fil. 
3.2 多 基因 累积 分 与 母亲 消极 教养 对 青少年 攻击 的 影响 

为 揭示 多 基因 与 母 灯 消极 教养 对 青少年 攻击 的 影响 ， 分 别 以 T1 青少年 攻击 行为 、T2 
青少年 攻击 行为 为 因 变 量 ， 性 别 为 控制 变量 ， 多 基因 累积 分 、 母 亲 消 极 教养 、 多 基因 累积 分 
x 母 亲 消 极 教养 交互 项 、 多 基因 累积 分 二 次 项 、 其 与 母 业 消极 教养 交互 项 为 预测 变量 进行 分 
层 回 归 分 析 。 

TL 青少年 攻击 行为 作为 因 变 量 的 分 析 结 果 如 表 2 所 示 。 母 亲 消 极 教养 显著 正 向 预测 Tl 
攻击 行为 ， 多 基因 累积 分 、 其 与 母亲 消极 教养 交互 项 对 T 攻击 行为 的 预测 作用 不 显著 ， 多 
基因 累积 分 二 次 项 对 Tl 攻击 行为 的 预测 作用 不 显著 ， 多 基因 累积 分 二 次 项 与 消极 教养 交互 
预测 TI 攻击 行为 。 简 单 斜率 分 析 ( 见 图 1) 显示 ， 母 亲 消 极 教养 水 平 较 高 (+1 SD) 时 ， 多 
基因 累积 分 二 次 项 对 Tl 攻击 行为 预测 显著 且 为 正 〈 二 次 项 : b= 0.05, t= 2.75, p= 0.006; 一 
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3 独立 样本 :检验 中 , 攻击 行为 .母亲 消极 教养 的 方差 不 齐 性 .为 校正 统计 偏差 , 采用 校正 后 的 自由 度 计算 t 值 CTabachnick & Fidell, 
2019) 。 此 处 报告 度 是 校正 后 的 自由 度 。 


次 项 :0 = -0.06, t=-1.92,p=0.055) ， 多 基因 累积 分 与 Tl 攻击 行为 间 呈 U 型 关系 ， 多 基因 
累积 分 高 或 低 的 青少年 均 会 表现 出 较 高 的 攻击 行为 ; 在 母亲 消极 教养 水 平 较 低 (-1 SD) 时 ， 
多 基因 累积 分 二 次 项 对 Tl 攻击 行为 预测 作用 不 显著 《二 次 项 : b= -0.03, t= -1.33, p = 0.18; 
— RI: b = 0.01, t= 0.17,Pp=0.86) 。 

#2 多 基因 累积 分 与 母亲 消极 教养 对 TI 攻击 行为 的 影响 


第 一 层 第 二 层 第 三 层 第 四 层 第 五 层 
AR? B p B p B p B p B p 
性 别 0.05”™ 0.22 <0.001 0.21 <0.001 0.21 <0.001 0.22 <0.001 0.21 <0.001 
NP 0.01™ 0.10 0.002 0.10 0.002 0.09 0.003 0.03 0.38 
MGPS —0.03 0.36  —0.03 0.37 = 0.03 0.31 —0.04 0.21 
MGPSxNP 0.00 —0.01 0.70 -0.02 0.60  —0.05 0.15 
MGPS? 0.001 0.04 = 0.22 0.03 0.42 
MGPS?xNP 0.007" 0.11 0.006 


TE: NP 代表 母亲 消极 教养 行为 ，MGPS 代表 多 基因 累积 分 ， 显 著 与 边缘 显著 的 结果 均 加 粗 显示 ， 下 同 。 
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图 1 多 基因 累积 分 与 母亲 消极 教养 对 T1 攻击 行为 的 影响 


T2 青少年 攻击 行为 作为 因 变量 的 分 析 结 果 如 表 3 所 示 。 母 亲 消 极 教养 显著 正 向 预测 T2 
攻击 行为 ， 多 基因 累积 分 、 其 与 母亲 消极 教养 交互 项 对 T2 攻击 行为 的 预测 作用 不 显著 。 多 
基因 累积 分 二 次 项 对 T2 攻击 行为 的 预测 作用 显著 〈 二 次 项 := 0.03, t= 1.98, p = 0.048; 一 
次 项 : = -0.02, t= —0.78, p = 0.44) ， 多 基因 累积 分 过 高 或 过 低 的 青少年 均 表 现 出 较 高 的 攻 
二 行为 〈 见 图 2) 。 多 基因 累积 分 二 次 项 与 消极 教养 显著 交互 预测 T2 攻击 行为 。 简 单 斜率 
分 析 显 示 ， 母 亲 消 极 教 养 水 平 较 高 G1 SD) 时 ， 多 基因 累积 分 二 次 项 对 T2 攻击 预测 显著 
HNE ZKH: b= 0.07, t= 3.20, p = 0.001; — RT: b =-—0.04, t=-—1.33,p=0.18) ， 多 基因 
累积 分 与 T2 攻击 行为 间 呈 U 型 关系 , 多 基因 累积 分 过 高 或 过 低 的 青少年 均 会 表现 出 较 高 的 
攻击 行为 ( 见 图 3) ; 母亲 消极 教养 水 平 较 低 (-1 SD) 时 ， 多 基因 累积 分 二 次 项 对 T2 Bee 
行为 的 预测 作用 不 显著 (二 次 项 : b = -0.02, t= -0.70, p = 0.49; 一 次 项 :b= 一 0.002, t = -0.06, p 


cu 


oT 


=0.95) . 4° 
表 3 多 基因 累积 分 与 母亲 消极 教养 对 T2 攻击 行为 的 影响 


第 一 层 第 二 层 第 三 层 第 四 层 第 五 层 
AR? B p B p B p B p B p 
性 别 0.04™ 0.21 <0.001 0.20 <0.001 0.19 <0.001 0.20 <0.001 0.19 <0.001 
NP 0.01™ 0.10 0.003 0.09 0.003 0.09 0.004 0.03 0.38 
MGPS 一 0.02 0.57 一 0.02 0.56 一 0.02 0.44 —0.03 0.31 
MGPSxNP 0.00 0.01 0.81 <0.001 0.99 —0.03 0.39 
MGPS? 0.004* 0.06 0.048 0.05 0.12 
MGPS2xNP 0.006 0.11 0.01 
0.1 0.1 
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图 2 多 基因 累积 分 对 T2 攻击 行为 的 影响 图 3 多 基因 累积 分 与 母亲 消极 教养 对 T2 攻击 行为 的 影响 


3.3 敏感 性 分 析 

起 对 任 一 单 基 因 、 母 亲 消 极 教养 对 TI1、T2 攻击 行为 影响 的 回归 分 析 结 果 显 示 ， 任 一 单 基 
.一 因 二 次 项 (AR28 < 0.004; |b|s <0.03, ps > 0.16) 、 其 与 母亲 消极 教养 交互 项 均 不 显著 (ARs < 
0.002; |b|s < 0.04, ps > 0.13， 见 附录 表 S1 和 表 S2) ， 表 明 多 基因 累积 分 在 提高 遗传 解释 率 和 
显著 性 方面 比 单 基因 研究 更 具 优 势 。k fold 交叉 验证 结果 显示 ， 模 型 预测 的 攻击 分 数 与 实际 
攻击 分 数 显 著 相 关 (rs>0.21,ps<1.46 x 10"), AES 10 次 后 > 值 仍然 稳定 显著 (rs = 0.20 
~ 0.24, ps = 3.91 x 10° ~ 2.87 x 1010， 见 附录 “2 k fold 交叉 验证 ”) 。 将 总 样本 随机 分 成 2 
个 子 样本 进行 的 内 部 验证 结果 显示 , 在 2 个 子 样本 中 , 多 基因 累积 分 二 次 项 与 母亲 消极 教养 
的 交互 作用 显著 或 呈 显 著 趋 势 〈 见 附录 表 S3 和 表 S4) 。 综 上 ， 本 研究 的 主要 研究 结果 获得 


验证 。 


4 本 研究 中 攻击 行为 得 分 的 分 布 是 正 偏 态 。 参 照 既 有 研究 (Harden etal., 2015) ， 对 攻击 行为 得 分 进行 对 数 转 换 并 重复 校正 前 的 
分 析 。 结 果 显 示 ， 各 统计 指标 显著 性 未 发 生变 化 ， 多 基因 累积 分 、 母 亲 消 极 教养 及 其 交互 项 对 T1、T2 攻击 行为 的 预测 作用 均 不 
显著 ;多 基因 累积 分 二 次 项 对 Tl 攻击 行为 的 预测 作用 不 显著 ， 对 T2 攻击 行为 的 预测 作用 显著 ， 多 基因 累积 分 二 次 项 与 母亲 消 
极 教养 显著 交互 预测 T1、T2 攻击 行为 。 需 要 注意 的 是 ， 对 攻击 数据 进行 转换 使 得 数据 弄 含 的 原始 信息 发 生 了 改变 ， 由 此 得 到 的 
可 归 分 析 结 果 中 参数 的 意义 解释 和 转换 前 有 所 改变 ， 反 而 不 利于 阐明 变量 间 的 关系 (Tabachnick & Fidell, 2019) 。 基 于 此 ， 我 们 
仍然 选择 在 正文 中 报告 不 校正 的 结果 。 
5 本 研究 选取 两 个 时 间 的 攻击 行为 进行 分 析 ， 可 能 涉及 发 展 变化 问题 ， 为 此 ， 我 们 控制 TL 攻击 行为 进行 了 探索 性 分 析 。 结 果 发 
现 ， 控 制 TL 攻击 行为 后 ， 多 累积 分 、 母 亲 消 极 教养 及 其 交互 项 、 多 基因 累积 分 二 次 项 及 其 与 母亲 消极 教养 的 交互 项 均 不 显 
著 ， 表 明 多 基因 累积 分 与 母亲 消极 教养 交互 预测 TI 和 T2 时 间 点 的 攻击 行为 水 平 ， 而 不 是 T1 到 T2 攻击 行为 的 变化 。 


4 讨论 
本 研究 采用 多 基因 累积 分 研究 范式 考察 了 多 巴 胺 系统 多 基因 与 青少年 期 两 个 时 间 点 的 
攻击 行为 间 的 关系 以 及 母亲 消极 教养 在 其 中 的 调节 作用 。 我 们 发 现 , 在 较 高 母亲 消极 教养 条 


件 下 多 基因 累积 分 与 青少年 攻击 行为 间 呈 U 型 关系 ， 多 基因 累积 分 过 高 或 过 低 的 青少年 均 
会 表现 出 较 高 水 平 攻击 行为 , 在 经 历 低 水 平 母亲 消极 教养 时 多 基因 累积 分 与 至 少年 攻击 行为 
的 关联 不 显著 。 


4.1 母亲 消极 教养 调节 多 巴 胺 系统 基因 与 攻击 的 关联 
多 巴 胺 系统 基因 与 攻击 间 的 U 型 关系 在 母亲 消极 教养 水 平 高 时 更 突出 ， 高 母亲 消极 教 
养 条件 下 ， 多 巴 胺 系统 多 基因 累积 分 二 次 项 对 T1、T2 攻击 行为 预测 均 显 著 ， 携 带 较 多 或 较 
少 与 低 多 巴 胺 活性 相关 等 位 基因 的 青少年 表现 出 较 高 的 攻击 行为 , 母亲 消极 教养 低 的 环境 中 
低 多 巴 胺 活性 相关 等 位 基因 与 个 体 的 攻击 行为 无 显著 关联 。 这 一 结果 与 GxE 交互 作用 理论 
和 相关 实证 研究 相符 (Beauchaine et al., 2007; Shanahan & Hofer, 2005; Tomlinson et al., 2022; 
Zhang et al., 2016) ， 表 明 高 消极 教养 环境 能 够 触发 或 者 加 剧 多 巴 胺 系统 多 基因 与 攻击 行为 
间 关 联 , 而 相对 积极 环境 (本 研究 中 即 低 母 亲 消极 教养 环境 ) 会 阻止 或 减弱 遗传 素质 的 表达 。 
7 母亲 消极 教养 促使 多 巴 胺 系统 多 基因 与 攻击 行为 间 关 联 增强 的 机 制 可 能 表现 为 以 下 过 
= 程 。 母 杀 消 极 教 养 作为 不 利 环境 影响 脑 内 基线 多 巴 胺 水 平 ， 促 使 多 巴 胺 系统 产生 诸多 变化 ， 
如 DRD2 受 体 密度 降低 、 中 脑 一 皮质 多 巴 胺 通路 过 度 激活 (Tielbeek et al., 2018) 。 由 于 青 
少年 期 多 巴 胺 能 敏感 性 以 及 对 压力 反应 性 增强 , 此 时 期 母亲 消极 教养 这 一 不 利 环境 因素 对 大 
脑 多 巴 胺 系统 的 影响 更 为 突出 。 因 此 ,多巴胺 系统 多 基因 累积 分 过 高 或 过 低 的 青少年 ,在 经 
历 高 水 平 母亲 消极 教养 时 , PFC 等 脑 区 功能 受到 基因 与 环境 的 双重 影响 而 受 损 ,导致 攻击 行 
为 风险 增加 。 同 时 ， 在 行为 层面 ， 母 亲 消极 教养 使 得 青少年 产生 对 攻击 或 敌意 性 的 模仿 、 愤 
怒 或 敌意 性 情绪 情感 增强 以 及 情绪 调控 能 力 受 损 等 (Gershoff, 2002; Goulter et al., 2020) , 
(S) 这 些 行为 方面 的 影响 与 过 高 或 过 低 的 多 巴 胺 功能 共同 作用 ， 导 致 青少年 攻击 行为 水 平 增高 

(Zhang et al., 2016) 。 

多 巴 胺 系统 多 基因 累积 分 与 攻击 行为 关系 符合 U 型 ， 与 低 多 巴 胺 活性 相关 联 的 多 基 
累积 分 过 高 或 过 低 的 青少年 均 表现 出 较 高 的 攻击 行为 ,这 与 已 有 神经 生化 研究 发 现 一 致 (如 
a ohne a 
和 过 低 的 多 巴 胺 含量 均 不 能 有 效 维持 皮层 兴奋 与 抑制 功能 的 平衡 ， 使 PFC 功能 减弱 ， 与 其 
他 脑 区 的 功能 连接 、 中 脑 一 皮质 多 巴 胺 通路 与 中 脑 一 边缘 多 巴 胺 通路 产生 异常 〈Arnsten， 
2009; Bertolino et al., 2009; Robbins & Arnsten, 2009; Seamans & Yang, 2004; Tian et al., 
2013) ， 从 而 限制 个 体 执行 控制 功能 、 情 绪 控 制 能 力 , 对 奖惩 的 敏感 性 增加 (Romeo, 2017) ， 
使 个 体 出 现 攻 击 行为 (Achterberg et al., 2020; Chester & DeWall, 2016; Padmanabhan & Luna, 


= 


加 


` 


2014; Raine, 2008; Rosell & Siever, 2015; Sadeh et al., 2019; Tielbeek et al., 2018) 。 尺 管 这 一 
研究 结果 需要 未 来 神经 生化 与 脑 科 学 领域 相关 研究 的 进一步 验证 , 但 是 该 结果 提示 ， 要 探究 


多 巴 胺 系统 基因 如 何 影响 攻击 行 
而 是 需要 关注 基因 本 身 的 复杂 性 
4.2 本 研究 的 创新 、 不 足 与 展望 

本 研究 的 重要 创新 在 于 采用 多 基因 累积 法 , 通过 累加 参与 多 巴 胺 代谢 、 传导 和 转运 的 多 
个 基因 位 点 , 揭示 了 多 巴 胺 系统 基因 与 青少年 攻击 行为 间 的 非 线 性 关系 。 这 对 于 理解 攻击 行 
为 的 复杂 生理 机 制 具有 重要 启示 。 尽 管 基 因 一 环境 交互 作用 仅 解 释 了 不 到 1% 的 变异 ， 但 对 
于 发 展 行为 遗传 学 研究 仍 具 有 重要 意义 (Evans, 1985; Hasan & Afzal, 2019) 。 基 于 基因 一 环 
境 交 互 作用 , 研究 者 可 关注 受 多 巴 胺 系统 基因 调控 的 生理 机 制 或 环境 敏感 性 特征 , 进一步 揭 
示 攻 击 行为 产生 的 生理 一 环境 机 制 ; 在 健康 决策 和 医疗 干预 中 , 基于 个 体 的 候选 基因 信息 或 
生理 机 制 特征 、 其 环境 敏感 性 ， 关 注 不 同 遗 传 素质 个 体 对 环境 (干预) 的 反应 性 差异 ， 从 而 
形成 更 具 针对 性 的 治疗 或 干预 方案 。 譬如, 为 具有 攻击 行为 且 对 环境 更 敏感 的 青少年 提供 积 


为 , 不 能 仅 依据 高 低 多 巴 胺 活性 判断 基因 “好” 还 是 “ 坏 ”， 
,考虑 多 巴 胺 系统 多 基因 组 合 与 攻击 行为 可 能 的 U 型 关系 。 


m 极 教养 环境 ， 能 够 在 一 定 程度 上 阻止 遗传 素质 的 表达 ， 减 少 攻击 行为 。 

ir 本 研究 揭示 青少年 期 多 巴 胺 系统 多 基因 对 攻击 行为 的 影响 表现 为 非 线性 模式 , 与 以 童年 
期 儿童 为 被 试 的 研究 发 现 不 同 (Davies et al., 2015; Thibodeau et al., 2015; Vrantsidis et al., 
= 2021) 。 这 与 大 脑 多 巴 胺 系统 的 发 展 动态 性 有 关 (Naneix et al., 2012; Tielbeek et al., 2018; Toth, 


2022) ， 也 与 青少年 期 个 体 需 进行 复杂 决策 、 学 习 抑 制 不 合理 行为 和 冲动 〈Tielbeek et al., 
2018), 且 该 时 期 大 脑 关 键 区 域 ( 如 , PFC) 尚未 发 育 完 全 , PFC 调控 功能 超 负 蓓 等 有 关 ( Toth, 
2022) 。 这 一 结果 提示 ， 青 少年 期 可 能 是 基因 一 环境 交互 作用 模式 发 生变 化 的 关键 期 。 未 来 
N 需 采 用 时 间 跨 度 更 长 的 追踪 设计 ,更 充分 地 探究 、 比 较 不 同年 龄 阶段 多 巴 胺 系统 多 基因 与 攻 
击 行为 间 的 关系 模式 ， 以 进一步 揭示 攻击 遗传 基因 机 制 的 发 展 动态 性 与 基因 作用 的 敏感 期 。 
同时 ， 本 研究 中 ， 多 基因 累积 分 二 次 项 、 其 与 母亲 消极 教养 的 交互 项 预测 T1、T2 的 攻击 行 
= 为 水 平 ， 不 能 预测 TL 到 T2 攻击 行为 的 发 展 变化 。 但 是 ， 这 并 不 意味 着 基因 一 环境 交互 与 
FE 攻击 行为 的 发 展 无 关 。 未 来 研究 应 深入 考察 遗传 与 环境 对 攻击 行为 及 其 发 展 变化 的 影响 。 
本 研究 存在 一 些 局 限 。 第 一 ， 本 研究 基于 攻击 行为 的 神经 生理 机 制 和 相关 理论 , 选取 了 
参与 多 巴 胺 代谢 、 转 运 和 传导 的 三 个 功能 性 位 点 ， 涵 盖 多 种 影响 多 巴 胺 活性 的 生理 过 程 ， 避 
免 了 数据 驱动 和 过 度 拟 合 的 风险 。 但 是 ,多巴胺 系统 基因 数量 众多 ， 本 研究 未 能 将 所 有 多 巴 
胺 系统 基因 或 攻击 行为 候选 基因 涵盖 在 内 ， 因 此 解释 率 仍 受 限制 。 未 来 研究 仍 需 采用 多 基因 
研究 范式 , 考察 多 巴 胺 系统 其 他 候选 基因 或 多 巴 胺 系统 与 其 他 系统 基因 间 的 联合 效应 ， 以 更 
全 面 地 揭示 攻击 行为 的 多 基因 遗传 基础 。 第 二 ， 虽 然 本 研究 进行 了 多 种 内 部 重复 验证 ， 为 看 
完结 果 的 可 靠 性 提供 了 证 据 。 但 随机 分 组 后 样本 量 减 少 可 能 会 降低 统计 检验 力 ; 本 研究 未 有 
使 用 相同 测量 工具 及 相似 特征 被 试 的 外 部 样本 ， 无 法 进行 外 部 验证 。 未 来 尚 需 在 其 他 样本 ， 
尤其 是 更 大 样本 中 进行 外 部 验证 。 第 三 ， 本 研究 只 关注 了 城市 常态 群体 的 青少年 , 研究 结果 
可 能 不 适用 于 高 风险 青少年 群体 ,如 低 社会 经 济 地 位 或 临床 样本 。 由 于 高 风险 群体 中 母亲 消 
极 教养 水 平 可 能 更 高 、 存 在 攻击 之 外 的 共 病 现象 等 (Davies et al., 2019; Gard et al., 2022; 


= 


Stocker et al., 2017) ,基因 与 攻击 行为 的 关联 强度 与 模式 可 能 不 同 于 常态 青少年 群体 。 第 四 ， 
本 研究 仅 考察 了 母亲 消极 教养 , 未 涉及 父亲 教养 。 父亲 教 养 行为 对 青少年 问题 行为 具有 重要 
影响 (如 ，Cao & Rijlaarsdam, 2022) 。 未 来 研究 可 同时 关注 父亲 教养 和 母 杀 教养 行为 ， 以 
更 全 面 揭示 父母 教养 行为 对 青少年 攻击 行为 的 影响 作用 。 第 五 , 本 研究 仅 关注 青少年 的 基因 ， 
但 是 考虑 父母 基因 可 能 与 教养 行为 存在 关联 以 及 父母 与 孩子 间 的 代 际 传递 问题 , 未 来 有 必要 
也 测评 父母 的 基因 ， 以 更 全 面 探究 基因 与 环境 对 青少年 攻击 行为 的 影响 。 


本 研究 为 攻击 行为 的 多 基因 遗传 基础 以 及 多 巴 胺 系统 基因 与 青少年 攻击 行为 间 存 在 非 
线性 关系 提供 了 证 据 。 在 青少年 早 中 期 ， 多 巴 胺 系统 基因 与 攻击 行为 呈 非 线性 关系 ， 并 且 高 
消极 教养 环境 能 够 加 剧 或 放大 二 者 之 间 的 关系 。 这 一 结果 提示 ，, 攻击 行为 干预 需要 从 减少 消 
极 教养 、 增 加 积极 教养 入 手 。 
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250014, China) 
Abstract 
Dopaminergic genes have been frequently found to be associated with aggressive behavior, 
but the results are inconsistent. One reason for the inconsistencies is there might be the U-shaped 
relationship between dopaminergic genetic variants and aggressive behavior. More specifically, 
evidence has suggested an inverted U-shaped relationship between dopamine activity and 
prefrontal cortex (PFC) function (a critical region related to aggression), with both dopaminergic 
hypofunction and hyperfunction, were related to poor PFC function. It is possible that the 
relationship between dopaminergic gene and aggression approximates a U-shaped function. 
However, such U-shaped relationship is rarely investigated in previous studies. Moreover, several 
concerns have been raised about the ignoring the polygenic traits of aggressive behavior when 
conducting gene by environment interaction (GxE) research using single loci. Therefore, the 
present study aimed to examine the interaction between dopaminergic genetic variants and 
maternal negative parenting on adolescent aggressive behavior by adopting the approach of 


multilocus genetic profile score (MGPS). 


Participants were 1044 adolescents (mean age 13.32 + 0.49 years old at Time 1, 50.2% 
females) recruited from the community. The adolescents completed two assessments with an 
interval of one year. Saliva samples, mother-reported parenting data and data on peer-nominated 
aggressive behavior were collected. All measures showed good reliability. The MGPS was created 
by COMT 1s4680 polymorphisms, DRD2 1s1799978 polymorphisms and DATI 1827072 
polymorphisms. Genotyping in three dopaminergic genes were performed for each participant in 
real time with MassARRAY RT software version 3.0.0.4 and analyzed using the MassARRAY 
Typer software version 3.4 (Sequenom). To examine whether negative parenting moderates the 
effects of MGPS on adolescent aggressive behavior, hierarchical regression analyses were 
conducted. 

The results found that after controlling for gender, maternal negative parenting was a 
significant risk factor for adolescent aggressive behavior, with higher negative parenting related to 
more aggressive behavior. The main effect of the quadratic term of MGPS on adolescents’ 
aggressive behavior was significant at Time 2, indicating a U-shaped relationship between MGPS 
and adolescent aggressive behavior. Moreover, the quadratic term of MGPS significantly 
interacted with maternal negative parenting in predicting aggressive behavior at Time 1 and Time 


2, respectively. Specifically, there was a U-shaped relationship between MGPS and adolescent 


aggressive behavior, indicating that adolescents with higher and lower MGPS exhibited higher 
levels of aggressive behavior when experiencing higher levels of maternal negative parenting. No 
significant effect of MGPS on adolescent aggressive behavior when experiencing lower levels of 
maternal negative parenting existed. 

This study provides evidence for the molecular mechanisms of multilocus genetic profile 
scores and gene-environment interactions in adolescent aggressive behavior. 
Key words aggressive behavior, maternal negative parenting, dopamine, multilocus genetic 


profile score, U-shape relationship 


附件 系列 敏感 性 分 析 


1 单 基因 效应 检验 
对 任 一 单 基因 与 母 杀 消极 教养 对 Th T2 攻击 行为 的 影响 进行 回归 分 析 。 结 果 显 示 ， 母 
灯 消 极 教养 对 T1、T2 攻击 行为 的 预测 作用 显著 ， 任 一 单 基因 及 其 与 母亲 消极 教养 交互 项 、 


任 一 单 基 


因 二 次 项 及 其 与 母亲 消极 教养 交互 项 


WARES OH $S1、S2) ， 表 明 多 基因 累积 


分 在 提高 遗传 解释 率 和 显著 性 方面 比 单 基因 更 具 优势 。 
表 S1 单 基因 与 母亲 消极 教养 对 TI 攻击 行为 的 影响 
第 一 层 第 二 层 第 三 层 第 四 层 第 五 层 
ARB p B p B p B p B p 
COMT 
Val158Met 
性 别 0.0S 0.22 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 
NP 0.01" 0.09 0.003 0.09 0.003 0.09 0.003 0.07 0.13 
COMT 0.01 0.72 0.01 0.71 0.04 0.34 0.04 0.34 
COMTxNP 0.00 —0.01 0.82 -9.004 0.90 -0.04 0.46 
COMT? 0.001 —0.04 0.31 —0.05 0.27 
COMT?xNP 0.001 0.05 0.38 
DRD2 
rs1799978 
性 别 0.0S 0.22 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 0.21 + <0.001 
NP 0.01™ 0.09 0.003 0.09 0.003 0.09 0.003 0.09 0.05 
DRD2 —0.04 0.16 —0.04 0.16  —0.07 0.16 —0.07 0.16 
DRD2xNP 0.00 一 0.01 0.72 一 0.01 0.71 —0.02 0.73 
DRD2? 0.00 0.03 0.50 0.03 0.50 
DRD2?xNP 0.00 0.01 0.91 
DATI 
rs27072 
性 别 0.0S 0.22 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 0.21 = <0.001 
NP 0.01" 0.09 0.003 0.09 0.003 0.09 0.003 0.05 0.27 
DATI1 —0.02 0.50 —0.02 0.50 —0.02 0.64 —0.02 0.63 
DAT1xNP 0.00 一 0.01 0.81 一 0.01 0.81 0.02 0.56 
DATI? 0.00 0.001 0.97 0.002 0.96 
DAT12xNP 0.002 0.07 0.19 


n 


TE: NP 代表 母亲 消极 教养 行为 ， 显 著 与 边缘 显著 的 结果 均 加 粗 ， 下 同 。 


表 S2 单 基 因 与 母亲 消极 教养 对 T2 攻击 行为 的 影响 


第 一 层 第 二 层 第 三 层 第 四 层 第 五 层 

AR? B p B p B p B p B p 
COMT 
Val158Met 
性 别 0.04°* 0.21 <0.001 0.20 <0.001 0.19 <0.001 0.19 <0.001 0.19 <0.001 
NP 0.01” 0.09 0.004 0.09 0.007 0.09 0.007 0.07 0.12 
COMT 0.04 0.17 0.04 0.19 0.06 0.18 0.06 0.18 
COMTxNP 0.001 0.03 0.30 0.03 0.28 0.02 0.75 
COMT? 0.00 一 0.02 0.58 一 0.03 0.54 
COMT2xNP 0.00 0.03 0.58 
DRD2 
rs1799978 
性 别 0.04°* 0.21 <0.001 0.20 <0.001 0.19 <0.001 0.19 <0.001 0.19 <0.001 
NP 0.01” 0.09 0.003 0.09 0.003 0.09 0.003 0.08 0.08 
DRD2 -0.04 0.16  —0.04 0.16 —005 0.29  —0.05 0.28 
DRD2xNP 0.001 —0.03 0.31 —0.03 0.30  —0.05 0.32 
DRD2? 0.00 0.01 0.83 0.01 0.83 
DRD22xNP 0.00 0.02 0.69 
DATI 
rs27072 
性 别 0.04°* 0.21 <0.001 0.20 <0.001 0.20 <0.001 0.20 <0.001 0.20 <0.001 
NP 0.01” 0.09 0.003 0.09 0.003 0.10 0.002 0.04 0.33 
DATI -0.04 0.27  —0.04 0.27 0.02 0.60 0.02 0.61 


DATIxNP 0.00 


—0.002 0.94 —0.004 0.89 0.03 0.42 


DAT1? 0.004 0.08 0.054 0.08 0.051 
DAT12xNP 0.002 0.08 0.13 


2 k fold 交叉 验证 
为 了 进一步 验证 模型 的 延展 性 和 预测 性 ,采用 MATLAB R2013b 的 regress 和 crossvalind 
函数 ， 进 行 k fold 交叉 验证 (k= 10; Marcot & Hanea, 2021; 张 文 新 等 , 2021) ， 其 中 9 组 作 


为 训练 集 ，1 


组 作为 测试 集 ， 回 归 模 型 为 y= c+ Bi x MGPS + Bo x maternal negative parenting 


+ B3 x MGPS x maternal negative parenting + B4 x MGPS x MGPS + Bs x MGPS x MGPS x 


maternal nega 


tive parenting + Bo x gender + e。 结 果 显 示 ， 模 型 预测 的 攻击 分 数 与 实际 攻击 分 


数 显著 相关 (T1 攻击 :r= 0.24, p = 3.45 x 1014; T2 攻击 :r= 0.21,p= 1.46 x 101) 。 为 了 测 
试 结 果 的 稳定 性 ， 对 这 一 过 程 重 复 进 行 10 次 , > 值 仍然 稳定 显著 (T1 Wa: r= 0.22 ~ 0.24, p 


= 3.91 x 10 
3 内 部 验证 


将 样本 按照 50% 随 机 分 成 两 个 子 样本 Cn sexi = 500, n rex = 544) 进行 


进行 内 部 验 说 


~1.74 x 100;T2 Wak: r= 0.20 ~ 0.22, p = 9.81 x 101? ~ 2.87 x 1010) 。 


u 
= 


分析， 以 
E。 两 个 子 样本 主要 变量 的 差异 均 不 显著 (性 别 : 92 = 2.23, p= 0.14; 多 基因 累积 


分 : 2 = 3.20, p = 0.67; 母亲 消极 教养 :t= 0.49, p = 0.62; T1 KE: t=-0.24, p = 0.81; T2 Wit: t 
=-0.14, p=0.89) 。 如 表 S3, S4 所 示 ， 多 基因 累积 分 二 次 项 与 母亲 消极 教养 的 交互 项 对 T1 
攻击 行为 的 预测 在 子 样本 1 Cp = 0.03) 中 显著 ， 在 子 样本 2 (p = 0.059) PS BA; X 
T2 攻击 行为 的 预测 在 子 样本 1 (p= 0.06) 和 子 样本 2 (p= 0.056) 中 呈 显 著 趋势 。 
对 于 Tl 攻击 行为 ， 在 子 样本 1 中 ， 母 亲 消极 教养 水 平 高 (> 1.90 SD) 时 ， 多 基因 累积 
分 二 次 项 对 TI 攻击 行为 预测 显著 且 为 正 Cb = 0.09, t= 1.97, p = 0.05) ， 多 基因 累积 分 与 T1 
攻击 行为 间 星 U 型 关系 ， 多 基因 累积 分 过 高 或 过 低 的 青少年 均 会 表现 出 较 高 攻击 行为 ;在 
经 历 低 水 平 母亲 消极 教养 (-1 SD) 时 ， 多 基因 累积 分 二 次 项 对 T 攻击 行为 无 显著 影响 Cb 
= 一 0.06, t=-1.72,p=0.09) ; 在 子 样本 2 中 ， 和 母亲 消 极 教养 水 平 较 高 (+1 SD) 时， 多 基因 
累积 分 二 次 项 对 Tl 攻击 行为 预测 显著 且 为 正 (b= 0.07, 1= 2.40,p = 0.02) ， 多 基因 累积 分 与 
TI BG {TAMU 型 关系 ， 多 基因 累积 分 过 高 或 过 低 的 青少年 均 会 表现 出 较 高 攻击 行为 ; 
母亲 消极 教养 水 平 较 低 C1 SD) 时 ， 多 基因 累积 分 二 次 项 对 Tl 攻击 行为 无 显著 影响 (b = 
-0.01, 1= -0.38,p=0.70) 。 

对 于 T2 攻击 行为 ， 母 亲 消 极 教养 水 平 较 高 (+1 SD) 时 ， 多 基因 累积 分 二 次 项 对 T2 攻 
击 行为 预测 显著 且 为 正 ( 子 样本 1:5= 0.059, t= 1.92, p= 0.055; 子 样本 2: b = 0.08, t = 2.66, p = 
0.008) ， 多 基因 累积 分 与 T2 攻击 行为 间 呈 U 型 关系 ， 多 基因 累积 分 过 高 或 过 低 的 青少年 
均 表 现 出 较 高 攻击 行为 ， 母 亲 消 极 教养 水 平 较 低 (-1 SD) 时， 多 基因 累积 分 二 次 项 对 T2 
攻击 行为 无 显著 影响 ( 子 样本 1: b =-0.04, t =-1.00, p = 0.32; 子 样本 2: b= -0.006,7= -0.18， 
p=0.86) 。 


qj 


表 S3 多 基因 累积 分 与 母亲 消极 教养 对 Tl 攻击 行为 的 影响 
第 一 层 第 二 层 第 三 层 第 四 层 第 五 层 
AR è B pp B p B p B p B p 95% CI? 

子 样本 1 

Hf 0.05°* 0.22 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 0.21 <0.001 [0.14,0.35] 
NP 0.004 0.06 0.19 0.06 0.20 0.06 020 -0.01 0.81 [008,0.06] 
MGPS 0.001 0.98 0.001 0.99 0.001 0.99 -002 0.75  [-0.06,0.05] 
MGPSxNP 0.00 0.01 0.79 0.01 079 -001 0.84 — [-0.07, 0.05] 
MGPS? 0.00 0.002 0.96 0.02 0.66 [-0.05,0.03] 
MGPS2xNP 0.01" 0.13 0.03 [0.005,0.10] 
子 样本 2 

性 别 0.05°" 0.23 <0.001 0.21 <0.001 0.21 <0.001 022 <0.001 0.21 <0.001 [0.18,0.40] 
NP 0.02" 0.12 0.005 0.12 0.004 0.12 0.005 0.06 0.20 Ss 0.03, 0.11] 
MGPS 0.05 0.24 0.05 0.26 0.06 018 —0.06 0.15 [-0.10,0.02] 
MGPSxNP 0.001 0.03 047 0.04 040 -007 0.17 [0.11,0.01] 
MGPS? 0.004 0.07 0.13 0.06 0.19  [—0.01, 0.07] 
MGPS2xNP 0.006* 0.11 0.059 [-0.001, 0.08] 

TE: MGPS 代表 多 基因 累加 得 分 ， 下 同 。 
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累积 分 与 母 杀 消极 教养 对 T2 攻击 行为 的 影响 
多 第 二 层 第 三 层 第 四 层 第 五 层 
AR? p B B B 95% CT 

子 样本 1 

生 别 0.04°** <0.001 0.18  <0.001 0.18  <0.001 0.18 <0.001 [0.11, 0.34] 
NP 0.005 0.14 0.07 0.16 0.06 0.17 0.001 0.99  [-0.07, 0.07] 
MGPS 0.98 0.001 0.98  -0.003 0.95 -0.02 0.73 — [-+0.07,0.05] 
MGPSxNP 0.001 0.04 0.41 0.03 0.50 0.01 0.79 — [+0.06, 0.07] 
MGPS2 0.002 0.04 0.39 0.02 0.67  [—0.04, 0.06] 
MGPS?xNP “0.007+ 0.11 0.06 —_[-0.002, 0.10] 
子 样本 2 

性 别 0.05*** <0.001 0.21 <0.001 0.22 <0.001 0.21 <0.001 [0.18, 0.41] 
NP 0.01° 0.008 0.12 0.007 011 0.01 0.05 0.31 [-+0.03,0.11] 
MGPS 0.45  —0.03 0.47 — —0.04 0.33 —0.05 0.28 [—0.09, 0.03] 
MGPSxNP 0.00 —0.02 0.70  —0.02 0.59 0.06 0.20  [-0.10, 0.02] 
MGPS? 0.006* 0.08 0.064 0.07 0.10 [—0.01, 0.08] 
MGPS?xNP 0.007* 0.11 0.056  [-0.001, 0.08] 


